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Abstract 

, The 0{N) Heisenberg and spherical models with interaction given by the long range hier- 

archical Laplacean are investigated. Two classical results are adapted. The Kac-Thompson 
solution [KTJ of the spherical model, which holds for spacially homogeneous interaction, is 
firstly extended to hierarchical model whose interaction fails to be translation invariant. Then, 
K*- I the convergence proof of 0{N) Heisenberg to the spherical model by Kunz and Zumbach [KZj 

■ is extended to the long range hierarchical interaction. We also examine whether these results 

. can be carried over as the size of the hierarchical block L goes to 1. These extensions are 

CN I considered a preliminary study prior the investigation of the model by renormalization group 

Tlj" ■ given in |MCGj where central limit theorems for the spherical (iV = cxd) model on the local 

. potential approximation (L J, 1) are then established from an explicit solution of the associate 

I nonlinear first order partial differential equation. 



1 Introduction 

Motivations. In the present work we initiate a geometric study of partial differential equations 
related to the renormalization group transformation (RGT) of a (i-dimensional A^-component hi- 
erarchical spin system in the limit as the block size L'^ goes to 1. 

For one-component hierarchical spin system, the evolution equation corresponding to L | 1 
limit, known as the local potential approximation, began to be investigated by Felder [F] who has 
constructed, in that seminal work, a family of global stationary solutions - nontrivial fixed point 
analogs of the corresponding (i-dimensional models with d > 2 and L'^ an integer > 1. Except by a 
preliminary study and numerical simulation of the corresponding A^-component spin equation by 
Zumbach jZj and a refinement of Felder 's approach by Lima [Lj, no much attention has been paid 
to the evolution equation model. 
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Very recently, the renormalization group trajectories of d = 4 dimensional hierarchical 0{N) 
spin model with block size L*^ = 2 has been investigated by Watanabe [W] (for previous investiga- 
tions, see \GK\ [K] and references therein). Starting from the uniform "a priori" measure supported 
in the iV-dimensional sphere of radius \/N, the critical trajectory has shown to converge to the 
Gaussian fixed point for sufficiently large A^. To control such trajectory, which starts far away 
from the fixed point, the exactly solved 0{oo) trajectory has been used together with two key 
ingredients: refiection positivity and the Lee- Yang property of single-site "a priori" measures. The 
former ingredient gives uniform convergence of 0{N) trajectories to 0{oo) trajectories. The latter 
property has been previously employed by Kozitsky [K] to establish central limit theorems for the 
hierarchical 0{N) spin model with block size L'^ > 2 and d > 4 (in his notation L'^ = 5 and 
2/d = A G (0, 1/2)) at the critical inverse temperature jSc and below. Watanabe's analysis, based 
in his joint work with Hara and Hatttori |HHW] on the critical trajectory for the hierarchical Ising 
model (A^ = 1), in contradistinction, deals with the borderline d = A case and does not require 
closeness to the fixed point. 

Although the analysis of the RGT with L'^ > 2 fixed is expected to be simplified considerably 
in the L [ 1 limit, none of the results on the critical trajectories can be carried to the limit as the 
above mentioned ingredients do not hold if L*^ is not an integer. To establish weak convergence of 
the hierarchical 0{N) Heisenberg equilibrium measure to the corresponding spherical equilibrium 
measure as A^ — > cxd in the local potential approximation (L | 1) an entirely new method of analysis 
has to be developed from scratch. 

The present investigation establishes some classical results on the 0{N) Heisenberg and spherical 
model with short range discrete Laplacean interaction replaced by the long range hierarchical 
Laplacean and examines whether they can be carried over as L goes to 1. Kac-Thompson solution 
of the spherical model |KTj . which holds for spatially homogeneous interaction, is presented in 
Section [2] and extended in Section [3] to hierarchical model whose interaction fails to be translation 
invariant. Kac-Thompson's asymptotic analysis is then applied to the moment generating function 
of the block spin random variable and normal fiuctuations are established for the spherical model for 
P < Pc- To our knowledge, such application and consequences are new. In Section IHthe convergence 
proof of 0{N) Heisenberg to the spherical model by Kunz and Zumbach [KZj is extended to the 
long range hierarchical interaction. The prove holds for the free energy and the moments generating 
function. These extensions are considered a preliminary study prior the investigation of the model 
by renormalization group. In a subsequent work |MCG] we establish central limit theorems for 
the spherical model on the local potential approximation from an explicit solution of the associate 
nonlinear first order partial differential equation given by (11. 7p with N = oo. 

Viscosity Limit. The hierarchical Heisenberg model on a box A C Z'^ of size n = L'^^ is defined by 
an 0{N) invariant equilibrium measure 

dui'^^x) = ^exp |-i {x,Ax)\f[daif'^ (x,) (1.1) 

where ) denotes an element of the configuration space fi„ = M x ■ ■ ■ x , A = J®I 

the tensor product of the n x n coupling hierarchical matrix J (see (13.21) ) with the N x N identity 
matrix I and ctq (x) the a priori uniform measure on the A^-dimensional sphere \xf = jSN of radius 

The invariance of J under the block spin transformation (13. ip allows to describe the laws of 
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equilibrium studying the dynamics of a recursion relation 

in the space of single-site "a priori" measures, which in terms of their characteristic functions 

0f)(^)= / dai''\x) exp(z^-x) , (1.2) 



reads 

Nk " V 2 



0r(-) = ^expf^A)(0K(L-/^.)) (1.3) 



with 7 = d + 2 and initial condition 

00 [z) = —^———wj^^ [N 2) = <^'q '{\z\) 1.4 

Here, exp (tA) is the semi-group generated by the A^-dimensional Laplacean A = d'^/dzf + ■ ■ ■ + 
d'^/dzff, Nk is chosen so that 0fc(O) = 1 holds for all = 1, . . . , and Ja(x) is the Bessel function 
of order a. Note that 4>\^\z) = (p\!^\\z\) depends only on the norm \zf = z - z and (3 is the inverse 
temperature. 

We shall now explain our title. Since ipit, z) = exp (tA) ipoiz) = exp(— [/(t, z)) satisfies the heat 
equation with initial condition iplO, ■) = ipQ, U satisfies Ut — Af/ + |Vf/|^ = with [/(O, ■) = Uq = 
— hiipo and the sequence u"^\x), A; = 0, 1, . . ., defined (with = —Nx) by 

(^f)(v/=iV^)=exp(-iV«f)(x)) 

can be obtained by solving a nonlinear heat equatioiJl] 

Ut - ^xuxx -Ux + 2xul = (1.5) 

up to time t = 1/2 starting from the initial conditions m(0, x) = L''-u)^\[L~^ x): 

wf^x) =m(1/2,x) -m(1/2,0) (1.6) 

(the solution at a; = is subtracted to satisfy v^^^'*(0) = e"^"^ ''^"^ = 1). This is the starting point of 
Watanabe's investigation [Wj (see also [HHW] ) on Taylor coefficients v'^l-, / > 1, of the sequence 

of functions ^^^''(C) = ~ '^'k^\VN Q / N around C = 0. As already recognized by Watanabe, 
the parabolic equation fll.Sp becomes a first order hyperbolic equation when N —* oo and the 

expansion in powers of 1/A^ of a trajectory fo^' the A^-vectorial hierarchical model is a 

singular perturbation about the corresponding trajectory for the spherical hierarchical model. 

The local potential approximation replaces the exponent 1/2 in (11. 3p by (L — l)/2 and the 
interval of time evolved by (11. 5p tends to when L | 1. As a consequence, defining u^^\t,x) = 



^Subindex t and x refer to partial derivatives with respect to independent variables. 
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u\!^\x) for t = k\nL and taking the limit L I 1, k ^ oo with t fixed, the recursive initial value 
problem (11.51) and (11.61) for (^^i^'* j turn into a genuine initial value problem given by 

- - nf) + 2x (nf))^ + - du^^^ + ui^\t, 0) = (1.7) 

with u^^^{0,x) = — \'n.(pQ^\\^—Nx)/N = Mq^'*. Comparing to (11.51) . (I1.7p includes three extra 
terms, the last one ensures u^^\t, 0) = for all t > 0, corresponding to the operations of dilation, 
multiplication and normalization performed between two consecutive evolutions of (II. 5p . Note that 
the stationary solution u*{x) = —x of (II. 7p corresponds to the Gaussian fixed point of (11.30 . In the 
second work of our series [MCG] we give a geometric description of the trajectory \^u^°°\t, x),t > O}, 
in the viscosity limit = ooj^t and above the critical inverse temperature. Our third investigation 
will address the solution [u^'^\t,x),t > O} of (11.70 as a (singular) perturbation about the critical 
trajectory [u'^^\t,x),t > O}. 

Statement of Results. Equilibrium laws of the model are described by the distribution of block 
spin random variable = n~"'^^'^'^^ Sj=i limit as n — > cx) with a properly chosen 7. The 

characteristic function with respect to the equilibrium measure (II. ip of the block variable with 
7 = (i + 2 reads 

$f)(z) = yexp^^L-W)/2^^..^j^^W(^) 

exp {ix ■ z) da^^\x) = (f^^\\z\) . 

The equilibrium distribution )(x) = {X^ < x) converges weakly in the thermodynamic 
limit n = L'^^ — > 00 to a^^''{x) = u'^^^X'^ < x) if (p^^''{\z\) is continuous at origin and con- 
verges pointwise to a continuous (at origin) function (f^^\\z\). Hence, the equilibrium distri- 
bution z/(^)(XT < x) converges weakly to the equilibrium measure of the spherical model if 

limTv^oo [^^^\vN \z\) 1 = y9(°°^(|2;|) exist for every point z and coincides with the correspond- 
ing characteristic function of the latter model provided (y9*^°°''(|z|) is continuous at 2 = 0. These 
statements, which is independent of which order the limits n — » 00 and N ^ 00 are taken, are 
proven in Section H] for 7 = c?, L'^ > 2 integer and /3 different from the critical inverse temperature 
Pc = Pc{d, L) of the spherical model. 

The following result on the moment generating function holds for admissible coupling matrices, 
including hierarchical matrix. 

Theorem 1.1 The finite volume moment generating function of the Heisenberg model 

n N 



1=1 j=i 



"^l/N plays the role of viscosity since it is in front of the Laplacean as in the hydrodynamic equation of incom- 
pressible fluid. 
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with admissible reflection positive sequence of coupling matrices A, converges 

lim el^\(3,z) = e{f3,z) 

n,N—^c<y 



to the spherical model moment generating function 6(/3, z) (see 1^2. 34\ ) ) as n, N goes to infinity in 
any order, for (3 < (3c{A) given by Ii2.21\) and uniformly in compact intervals 0/2; G M. 



2 Spherical Model 

We review the solution and some basic properties of Berhn-Kac model with a positive definite 
coupling matrix A satisfying a condition stated in (12.121) . Formulas written in this section are inde- 
pendent on whether translational invariance holds and are, in addition, suitable to the hierarchical 
coupling matrix investigated in the next section. We shall make most of those expressions explicit 
by choosing A the usual discrete Laplacean, denoted here by —A. The same symbol will be used 
for the hierarchical Laplacean in Section [31 



2.1 The Free Energy 

Given (3 > and a positive coupling matrix J = [Jj^j^^^^, the spherical model of Berlin and Kac 
|BK] associated with f3 and J is defined by the partition function 

Qn{P,J) = ^ j dan{x]y/n) exp (a;, Ja3)| (2.1) 

n 

where {x, y) = ^^a^j Vi denotes the inner product in M", 



i=l 



danix;r) =5i\\x\\-r)Y[dx^ (2.2) 



i=l 
2 



the uniform measure on the sphere S„(r) = ja:; G M" : ||a3|| = {x, x) = r^} of radius r and 

f 27r"/^n(""^)/2 
Sn = J dan{x; y/n) = Y{n/2) ^'^'^^ 

is the surface of the sphere S„ (y/n). 

The most common choice of coupling matrix J is given by the discrete Laplacean, — AaEI on a 
d-dimensional hypercube A C Z"^ of size n = L'^ with periodic boundary condition: 

-(AaA= {h-f,) = m- E (2.4) 

j:\i-j\ = l j:\i-j\ = l 

with the summation over the lattice sites j which are at unit Euclidean distance from i G A. 



^For simplicity, we drop the subindex of Laplacean Aa if no confusion exists. 
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Remark 2.1 The Berlin-Kac model incorporates essential features of the ferromagnetic Ising 
model, exhibits a phase transition and has the advantage to he exactly solvable in any dimension. 
See IBK] for an extensive discussion on the thermodynamic properties above and below the critical 
temperature. The phase transition on the spherical model is of the same nature of that observed 
in the free Bose gas in which condensation of a single mode occurs ( see e.g. [P]). Disordered 
mean spherical model and equivalence of ensembles has been investigated by Pastur (see IKKPSf 
and references therein). See also Perez-Wreszinski-van Hemmen fPWH^ for disordered spherical 
models. 

To solve the spherical model it is convenient to introduce an auxiliary expression 



= Y j exp|-^(a;, (J-/i)x)| (2.5) 

which, after changing variables r = y/ns and x = sy, can be written as 

ds 

In = Vnl — exp{n /;,„(s)} (2.6) 



where 
and 



hn{s) = ^s^ + \ns + fn{s') (2.7) 



fn{s^) = ^\nQn{s\j) (2.8) 
is the finite volume free energy of the spherical model. 

In may be think as the grand-canonical partition function with the Lagrange multiplier yU < 
playing the role of a chemical potential. The function /„ can be integrated 

n("-i)/Vdet ( J - /i) 

and equations (12. 6p and (12.91) used to evaluate the free energy when n ^ oo. 

For instance, if J is given by (12. 4p . Fourier spectral analysis (see [D]) can be used in (12. 9 p 
together with Stirling's formula V{n/2) ~ (n/2e)"'^^ to write 

lim -ln4 = -^-^Eln(-A-^) (2.10) 



where 



Eln(-A-/i) = lim ^ V In (cu (27rm/L) - /i) 



-L/2<mi<L/2 



I d^k \n{u{k)- ^i) (2.ir 
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and 

d 



uj{k) = 4 sin 



2 

bill 

1=1 
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Note that, as n oo, —A is unitarily equivalent to an operator of multiplication by uj{k) in the 
space L2 ^[— vr, vr]'^ , of square integrable functions / : [— vr, vr]"^ — > C. 

Let us now state our assumption on J and explain the probabilistic notation E (■) in (12. 101) . 

Definition 2.2 A sequence A = {^n}„>i of coupling matrices (n indicates the order of An) is an 
admissible sequence if each An is nonnegative ({An)^j > positive definite real symmetric matrix 
and 

KfiA)= hm -TrfiAn) (2.12) 

n— >oo 11 

exists for every continuous bounded function f. 

We require, in addition, that 1 = {l/^/n) {1, . . . ,1) is an eigenvector of An with associate 
eigenvalue 0. 

From now on, only admissible sequences of coupling matrices will be considered. By definition, 
E/(A) = hm i^ / (Af)) = hm / t/p„(A) /(A) = / rfp(A) /(A) (2.13) 

i=l 

is the expectation with respect to the empirical distribution p which is the weak limit of the 
integrated density of eigenvalues A|^"'', . . . , An"'' (counting multiplicity) of An- 

1 " 



with X[a,fe)(A) = 1 if A G [a, 6) and X[a,fe)(A) = otherwise. 

The empirical measure p associated with —A is absolutely continuous with respect to the 
Lebesgue measure dp{X) = p'{\)d\, 

p\\) = l\im-s{i-A-\-^s),,') 

exists for almost every A and is supported in the interval [0,4(i]. For d = 1, we have explicitly 
p'(A) = 1 /(7rV4A-A2). 

On the other hand, Laplace's asymptotic method applied to (12. 6 p gives 



^^P{^h{s)}{l+0{l/n)) (2.15) 
with s the value at which h{s) = lim /i„(s) attains its maximum value. The existence of a unique 

n—*oo 

strictly positive maximum s follows from certain facts that are independent of J considered. Writing 
Qn{s'^) = (5n(s^, J), we have 
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1. Qn{0) = 1 

2- Q'nis^) = ^ I dcTnix; y^) (sj, Joj) exp | ^ (aj, Jaj) | 
3. 

< Q':,{s')Qn{s') - {Q'nis'f = ^J d^nix; v^) da^y; V^) 
X iix,Jx) - {y,Jy)Yexp { [{x, Jx) + {y,Jy)] 



If J satisfies (12.121) then || J|| < c is bounded (when J = —A, for instance, we have 

1 1 9 

< - (a;, Jx) < - \\J\\ \\x\\ = 2dn , 

in view of (12. 2p and ||J|| = sup^g[_^ co'(fc) = M) and together with the mean value theorem, we 
conclude fn^s^) = f'nis^) holds for some < s < s, with 

and 

uniformly in n. {fn{s^)} is a sequence of convex functions, uniformly bounded in every compact of 
and there is a subsequence l/n (s^)} such that lim (s^) = /(s^) exists with /(s^) convex 

and differentiable at almost every G M+ (See |KT] and Section 3.1 of |KKPSj ). As a consequence, 
the maximum of h is attained at 

for every fi < 0. 

Equating fl210|) and (12151) together with ([22D yields 

+ lns-+ /(.-2) = -i - iEln(J - ^) 

which, by the inverse function theorem, can be differentiated with respect to fi: 

f{fi)=E{J-i2y^ (2.17) 
in view of h'{s) = 0. When J = —A, the equation reads 

= 7:!— J d'^k — i . (2.18) 
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Solving (12.171) for /i = /i(s^), with replaced by an arbitrary positive number /3, substituting 
back to the previous equation gives the free energy of the spherical model 

= - In (v^) - IeIu (J - (2.19) 

(for J = — A the last term is given by (12.111) ). 



Remark 2.3 The sum-rule \2.11^ can he obtained directly from the grand-canonical partition func- 
tion ( f^) (see e.g. f^) 

2 2d\nln 
s = lim 

n^oo n OjJ, 

and this expresses the equivalence between different ensembles of this model. As a function of 
fi G (—00,0], Hn{n) = {2/n) dlnln/dfi is convex, monotone increasing with if„(— oo) = and 
lim^io -f^n(yu) = oo. If projects on the invariant subspace of —A associated with X = 0, it 
follows 

si = E^o(-A -/iJ)"^ 



> ^—r I d'^k ^ = s^-E(-A)-^ (2.20) 

The eigenvalue is said to condensate if > 0. According to the above inequality, Sq > provided 
P = s"^ > K (— A)~^. Since i2.11\) has a unique solution /i = /i(s) with Sq = for < E (— A)~"^, 
the spherical model exhibit a phase transition of Bose-Einstein type whenever E,{—A)~^ is finite, 
i. e. when d > 2. Note that, for k near 0, u{k) ~ \k\^ and 

F \ Inl/e ii d = 2 

has no limit for d < 2. For any sequence A of admissible coupling matrices, the critical inverse 
temperature of the spherical model is defined by 

l3c{A) = EA~^ (2.21) 



e 



2.2 Moments Generating Function 

We state the main result of this section. 
Theorem 2.4 The block spin random variable 

1 " 

X„ = ^Vxi. (2.22) 

of a spherical model with admissible sequence J = (Jn)>i of coupling matrices, converge in distri- 
bution to a Gaussian variable of zero mean and variance — l//i where fi = ^{(3) is the solution of 



( fO?P with r = 13. 
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Remark 2.5 Since fi approaches from below as f3 ^ {3^, the variance — diverges and Theorem 



2.4 does not hold anymore. Theorem 2.4 holds for cases (e.g. for no translational invariant J's) 



in which the general result of Newman ^Nj on Gibbs measure satisfying FKG property with finite 
susceptibility x cannot be applied . 



Let 



Qn{P,z,h,J) = ^ j dan{x;y/n) exp f^-^- {x, Jx) + z {h,x)^ (2.23) 

be the partition function of the spherical model including a magnetic field h = {hi, . . . , hn). Since 
Qn{P, z, h, J) = QniP, —z, h, J), fl2.23p is an even function of z and Qn{P, 0, h, J) = Qn{P, J) is 
given by (12. ip . One easily verifies that the ratio of partition functions 

with h in fl2.23p given by the ra-component unit vector 1 = {l/^/n) (1, . . . , 1), generates the mo- 
ments of Xn- 



r.p^ [ dan{x;y/^) XPexp!.-^(x,Jx)\ 

/ dan[x] y/n) exp < -- (x, Jx) 



and the same procedure of the previous subsection can be used to evaluate the moment generating 
function B„. Despite of fact that h in this case goes to when n ^ 00, the ratio (12.241) converges 
to a nontrivial (7^ 1) function of z as we shall see in the following. 

Repeating (12. 5p - (12. 9p with J„ replaced by an auxiliary function for (I2.23P yields 

If"- c_i 

Kn = — ]^c?Xi exp<—{x,{J~K)x) + {h,x) 

2^/2-1 Y(nl2) (If 1 

exp < - [ h, — -h 



n(--i)/Vdet ( J - /t) \2\'J-kI 
/"°° ds 

^/n I — exp{n 5(„(s)} (2.25) 
Jo s 



with 



where the free energy 



9n{s) = -S^ + Ins + fn{s^,s) 



fn iP, z) = -ln Qn {P, z; h, J) 
n 



is a smooth function of {P, z) such that (s^, 0) = /„ (s^) is given by (12.81) . 
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We continue from the last two equations of fl2.25p : 

1 

n^oo Tl 



lim-lnir„ = -i-^Eln(J-K) + ^E;,(J-K)~' 



+ Ins + /(s^s 



7^2 -\ 



(2.26) 



where s = s{k) is the solution of equation (12.161) with /i and /'(s^) replaced by k and the derivative 
of /(s^, s) with respect to (recall /(s^, s) is an even function of s): 



-1 



(2.27) 



and, if Ph = hh / \\h\\ denotes the projector in the h direction and {Xj, Ej}^.^-^ are the spectral 

n 

elements of J: J = \j Ej, we have 



||^||2 ^ 1 " 1 

Eh {J - = lim ^^Tr- Ph = lim - V \\Ejh\ 

n^co n J — K n^oo n ^-^^ A,- — K 

When J = — A, Fourier analysis gives 

Efc (-A - Kf^ = ^— r / d^k- 



(2.28) 



=1 



h{k) 



with 



{2nf J[-n,n]y '"uj{k) - K 
h{k) = hx exp {ik 



X) 



Since 1 = (l/i/n) (!,...,!) is the unique eigenvector of (12. 4p associated with the eigenvalue 0, only 
the zero-mode contributes to the expectation (12.281) with ^ = 1: 



Eh (-A - k)-^ = {h, (-A - k)-^ h) 



UK 



(2.29) 



Assuming 1 an eigenvector of J (orthogonal to the complementary space in view of J = J^) with 
associate eigenvalue A = the same result holds with —A replaced by J. We shall continue our 
calculation of the generating function 0„ with J satisfying the assumptions of admissible coupling 
matrices. 

The free energy function is obtained by equating the two lines of (I2.26p up to order 1/n and 
proceeding as in equations (I2.16I) - (I2.19I) : 



T) Z 

nfn (/3, z) = --{l3t^ + In (e/3) + E In ( J - /t)) - — + 0(1) 



where k„ = k„ (/3, z) is the solution of 



1 z' 
f3= + 



—A — K ' nn"^ 
11 



(2.30) 



(2.31) 



and the order 1 term is independent of z by (12.151) for Kn together with (12.321) below. 

Denoting by Wn{n,z) the r.h.s. of (I2.3ip . w„ is a decreasing function of k for n < with 
Wn{—oo,z) = and \imf^^QWn{n, z) = oo. By the imphcit function theorem, the solution k„ {P,z) 
of (]2.3ip is the unique real analytic function of z'^ in a neighborhood of z'^ = 0. Consequently, 
Kn {P, 0) = fi for every n, lim^^oo t^n (/?, z) = ^ uniformly in z, where is the solution of 
fl2Tn) and 

K = /i+ + ( — 1 (2.32) 
n V ^ / 

by Taylor theorem, where 

c= [E(l-J/^)-2]-' 

is obtained by plugging the resolvent equation 

1 111 
(K-/i) 



J — K J — J — K J — fl 

into (ICT]) together with fl2:T7D and (l2:32|) . 



Proof of Theorem 2.4 Substituting exp {nfn{f3, z)) in the numerator and denominator (with z = 0) 
of (12.241) . taking into account that 0(1) term in (I2.30p does not depend on z, gives 

e„(/3,.) = exp (-nf (K - ^) - ^Eln^ - il) (l + o (i)) (2.33) 

and it suffices to verify that 

In e(/3, z) = lim In e„(/3, z) (2.34) 

n— >oo 

converges uniformly in compacts of z and for every /5 < /?c to a limit proportional to z"^. The limit 
(I2.34P exists by the same reasons employed to show existence of lim^^oo /n(s^). Plugging (I2.32p 
into (1^:^ together fITTTp gives 

lne(/3,2) = ^ c/5-cE 



J-/i fij 2/i(/5) ■ 

□ 



3 Hierarchical Spherical Model 
3.1 Hierarchical Laplacean 

Paiva and Perez [PPj have investigated the semi-groups generated by d-dimensional Hierarchical 
Laplacean —A in the presence of disorder by using, for the first time, spectral analysis. Although 
—A has discrete spectrum, they have shown that exp (tA) 6o, with 6o localized at origin, diffuses. 
We quote |Kr] and references therein for spectral localization in hierarchical Anderson model. Here, 
in order to apply the limit theorem established in Subsection 12.21 we need the spectral theorem for 
homogeneous hierarchical matrices. We extend the work of Watanabe [W] to arbitrary spectral 
dimension d (see Remark [3.40 . 
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Given integer numbers L, K > 1 and d > 1, let 

Aj^ = {0,1,...,L^ - l}'^ C 

be a hypercube with cardinality \Ak\ = L'^^ = n. Let = {9i, . . . , 9k) denote the coordinates of a 
point i G Ak written in the L'^ base 



K 

I = 

k=l 



From now on we use these coordinates to index components of a vector u = (ue) in M^^. 
Let B : M/"^ — > M^^^'-i be the block operator 

6»G{0,lv,-f'-^-l}'* 

and let B* : M^^-i — > M^^ be its adjoint 

{B*v,u)^^ = {v,Bu)^^_^ 

with respect to the inner product {u, w)^^ = ue we'- 



{B%] 

Define a real symmetric matrix in M" 

K 



e 

1 



J = Y^L-^\B*f B^ . (3.2) 



k=l 



Proposition 3.1 The associate quadratic form of J for L'^ = 2 and d = 2/ (a — 1) gives the 
hierarchical energy 

K 2^-*= 

-H = J]2-"'=^(5,,,)' 

k=l r=l 

Sk,r = ^ 

(r-l)2'=<jr<r2'= 

introduced by Dyson [DJ in his study of the Ising model with 1/ \i — jl'^ -interaction. 
Remark 3.2 Note that a{d) = {d + 2) /d ranges from 2 to 1 as d varies from 2 to oo. 
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Proof. Setting L'^ = 2, we have Sk^r = 2^^^^ (^B^a)^ for some 6 G P^K-k and 

K K 

k=l k=l 

as claimed. 



□ 



We require the hierarchical Laplacean —A satisfies 

-A1 = (3.3) 
(Theorem 13.31 below shows that is also a simple eigenvalue). Writing 

J = A + /io/ 
together with (13. 3p . (13. 2p and (13.11) . we have 

/^o = (1,^1)a.- = ^flL--' (^'1'^'i)a._. = (3-4) 

k=l k=l 

and 

K 

- ^ = ^1-^^ [-{B*f B^ + I) (3.5) 
fc=i 

generates a stochastic semi-group. 
Now we observe that 

B'' {B*f = BB* = I (3.6) 

holds for every k = 1, . . . , K and 

Pk = {B*fB^ (3.7) 

is an orthogonal = P^ projection matrix P^ = P^. on the subspace of vectors in R^^ which 
assumes constant value over blocks of size L'^''. It follows from (13.61) 

PjPk = {B*y B^ {B*f B^ = {B*y B^-^B'' = Pj 

p^p^ = [B*f B^ {B*y B^ = {B* f {B*y-^ B^ = Pj (3.8) 
hold for any j > k and we have the following inclusions 

Pk < Pk-1 <---<Pi<Po = I (3.9) 

in the sense that A < B if, and only if, {u, Au)^^ < {u, Bu)^^ holds for all u G R^^. 
Let 

Qk = Pk-Pk+i (3.10) 

for k = 0,1, . . . , K — 1 and 

Qk = Pk 

be the block fluctuation operator. 



14 



Theorem 3.3 (Spectral) The collection {Qk}k=o of n x n real orthogonal projection matrices 

QjQk = SjkQk (3-11) 

are the spectral partition of unit 

K 



k=0 

and 

K 

/(-A) = 5^/(Afc) Qu (3.12) 

holds with 

A. = (3.13) 

for any continuous function f : [0, 1/(L^ — 1)] — > M. It follows that —A is a positive definite 
matrix where X^, k = 0, . . . , K — 1, is an eigenvalue of multiplicity L'^(^^''\l — L^^) and Ax = a 
simple eigenvalue. 

Proof. The prove is essentially given in fW]. By (13.101) and (13. 8p 

QjQk = {Pj — Pj+l) {Pk — Pk+l) 

= Pj {Pk — Pk+l) — Pj+l {Pk — Pk+l) 

= {Pj - P^) - (^.+1 - P,+i) = 
for any k < j < K and the same holds for j < k < K. For j < k = K, 

QjQx = {Pj — Pj+l) Pk = Pk — Pk = 

and for j = k 

QkQk = {Pk — Pk+l) {Pk — Pk+l) = Pk + Pk+l — 2Pa,.+i = Qk ■ 

By definition, 

K K-l 

Qk = {Pk — Pk+l) + Pk = Pq — Pk + Pk = I ■ 



k=0 k=0 

Finally, by ([33]), ([32D and fl310D . we have 



K 

-A = J2l''H-Pj + i) 

K j-1 

j=l k=0 
K-l / K \ 

= Y.[Y. L-^'jQk + O-QK (3.14) 

k=0 \j=k+l J 
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which gives (13.121) with /(x) = x. It follows by (13.111) that (13.121) holds for any polynomial and, by 
Weierstrass approximation theorem, for any uniformly continuous function. 

Since projects on vectors in which are constant over disjoint blocks of size L'^^ , the rank 
of Pk is 

rankPfc = L"^^^-^) . 

By definition (I3.10p together with the inclusions (13.91) . the rank of the block fluctuation projector 
Qk is 

rankQ, = L'^^-^^ - L^C^-'^-D (3.15) 

for /c = 1, . . . , i^' — 1 and 

rankQii- = 1 
and these concludes the prove of Theorem 13.31 

□ 



Remark 3.4 The spectral measure fi^ associated with the vector 60 = {^e,j)jiz\j^ defined by 

POO 

{5eJi-A)5g)= / f{x)d^ifix) 



for every bounded Borel function f : M — > C, is given by 

L'^ — 1 

ld{k+l) ^ ~ '^^ ' IdK 



d^^ {x) = j-dik+l) ^ ~ ^-t) + JdK^ (^) 

A:=0 



by inspection of the matrix elements Qeg', and is independent of 6. As n = L'^^ tends to infinity, 
is the unique weak-* limit point of the corresponding empirical measure ^2.14 ) (see Theorem 
1.2 of Kritchevski IKr^ ). The number 

m In t 

is called spectral dimension of —A. 



3.2 The Free Energy 

To compute the free energy (12.19^ of the spherical model associated with j3 > and hierarchical 
Laplacean matrix J = —A we need to evaluate the expectation E with respect to the empirical 
measure of eigenvalues of —A in both the last term of the r.h.s. of (12.191) and in the implicit 
equation (I2.17P for /i = /i(/3). 

Using Theorem 13.31 together with (12.131) and linearity of trace, we have 

1 1 ^ 

.Tr/(-A) = ^/(A,)TrQ, 



]^dK •' > l^dK 

k=0 



(l-L-^)X;^-* /(A,) + ^/(A^) 



fc=0 



16 



in view of the fact that the eigenvalues oi Qk are and 1 together with fl3.15p . Hence, the 
subsequence {fnKi.^'^)} Km ^"^^^ energy with riK = L'^^ converges to (12.191) where 

E In (- A - /i) = (1 - L-") J2 L-"' In ( - ^ (3. 16) 

fc=0 ^ ^ 

and fj, = fj,{P) solves 



including the 0-eigenvalue contribution po = ~ f^^) ^ which, by Remark I2.3[ may have 

macroscopic occupation. 

Analogously to (12.201) . we have 

oo 

> /?-(i-l-^)(l2-i)^l-('^-2). 



k=0 

(1 - L-^) (L2 - 1) 



1 - L-^+^ 

which is strictly positive provided d > 2 and (3 > f3c{d, L) where 

, , (1 - L-'') (L2 - 1) 

is the critical inverse temperature of the hierarchical spherical model. 

Remark 3.5 The geometric multiplicity (L^ — 1)L'^^^~''~^^ of each eigenvalue \k, A; = 0, . . . , K — 1, 
of the hierarchical Laplacean (SJ^ can be lift by fluctuation projectors Qk,e depending on the index 
9 G Ax-k, such that 

Qk,eQk,e' = ^ee'Qk,e , 
and a nonhomogeneous Laplacean can be defined by 

K 

—A^^ = Afc,6» Qk,e 

k=o e&KK-k 

with 

Afc,0 = c Afc exp {Xk,e} 

where Xk is given by l{3.13\) . {Xk^e} chosen according to a common probability distribution P with 
mean with = Eexp (Xk^e) < oo. In this case, we have 

^Trf{-A) = (l-^-^)i;^"*Z5(L) E /(^A.exp{X,.,,}) + ^/(A^) 

fc=0 O&Kx-k 



(1 - L-") 5^L-*E/(cAfcexp{X,,o}) = E/(-A) 



fc=0 
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for almost every {X^ o} with respect to distribution F, by the law of large numbers. Here E/ denotes 
the expectation with respect to the product measure dp{X)F (dx) with p the empirical distribution 
relative to the eigenvalues {Xk} O'^d P the common distribution of variables X^ q. Expressions 
^2.19\) . ^3.16\) and \3.1''/\) are obtained accordingly. For instance, 



3.3 Continuum Hierarchical Laplacean 

Hierarchical Laplaceans have discrete eigenvalues. We shall now consider a continuous version 
obtained by a limit procedure. 

We shall take L [1 simultaneously to K ^ oo maintaining KliaL fixed equal to C G M+ U {oo}. 
Equation (13.171) . for instance, reads 

K{L) / r _ 1 X -1 



d / ^ . — e-'^^dy 

io exp(-2y)-exp(-2C)-2/i 

(l-oxp(-2C))/2 . 

dp{X) (3.18) 

A — p 

where dp{X) = p'{X)dX is absolutely continuous with respect to the Lebesgue measure dX with 



if A G [0, (1 — exp (— 2C))/2] and otherwise. So, the empirical measure for the eigenvalues of 
the hierarchical Laplacean —A with L I 1 converges provided K = K{L) increases faster than 
C (lnL)~^. We take C = oo, for simphcity. 

Accordingly, equations (I2.19p . (13.161) and (I3.17P holds with replaced by A^ {L — 1) and 

po > /5 - lim ^ ^ ^ ' 



Lii (1 - L-'^+2) (L - 1) 
2d 



d-2 

has a strictly positive limit provided d > 2 and /? > Pc{d) where 



2d 

m = j-^ . (3.19) 



Note that 

■1/2 f"\,,. 2d 



El 



-A)-i = / A-Mp(A) = 2'^/^^ / xdl-i--idX 
Jo 2 Jq 



d-2 



18 



and a phase transition of Bose-Einstein condensation type occur at the critical inverse temperature 
/?c = E (— A)~^ as in the spherical model with the usual Laplacean interaction. 

We now compute integral (I3.18P with C = oo and (1 = 4, for comparison purposes. For /i < 0, 
we continue 

E (-A -,/)-= 8 f = 4 (2, in (l - i-) + l) 

and equation (I2.17P with A given by the hierarchical Laplacean at d = 4 and = /?, reads 

l_^._2;.l„(l-i-) . (3.20) 



4 Convergence to Spherical Model 



4.1 The 0{N) Heisenberg Model 

The partition function of the O {N) symmetric Heisenberg model in a ci-dimensional cubic box 
A C Z'^ of cardinality n = L*^ is given by 

A) = ^ f exp (-§ {x,Ax)\fldalf'^ (x,) (4.1) 

where x = {xi, . . . , x„) is a n-tuple with each Xj G M^, A = J ® I is the tensor product of a n x n 
coupling matrix J with the N x N identity matrix I and cTq^'' (dx) = (^dx; ^ Sn is the 

2 / 

"a priori" uniform probability measure on the A^-dimensional sphere \x\ = N of radius yN with 

n 

surface area Sn given by (12.31) . The inner product on M" ®]R^ is denoted here by {x, y) = Xi-yi. 

i=l 

The expected value with respect to the Heisenberg measure vi!^'^ (see (11.11) ) is defined by 

(^>4"' = ^ ^(-) {-f n^-o"' (-,) 

and from here on, A in (■) (n) = (■) (N)(p,A) is assumed to be a sequence of admissible coupling 
matrices in the following sense. 



Definition 4.1 A sequence A = of coupling matrices is an admissible reflection posi- 

tive sequence if each An is admissible in the sense of Definition \2.S\ and (■) {n){(3,A) satisfies the 
reflection positivity condition lli4-4\ )- 



Kunz and Zumbach |KZ] have devised a way of proving convergence of the free energy of the 
Heisenberg model to the spherical model for nearest neighbor interactions. We shall first show 
that their method holds for admissible reflection positive coupling matrices and prove Theorem li.il 
afterwards. 
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Theorem 4.2 The finite volume free energy of the Heisenberg model 



fi''\P) = ^^^Zif\(3,A) (4.2) 
with admissible reflection positive sequence of coupling matrices A, converges 

lim = (4.3) 

to the spherical model free energy Ii2.19\) as n, N goes to infinity in any order. 



The original proof by Kac and Thompson |KTj asserts that (14.31) holds for all coupling matrix J 
satisfying translation invariance Jij = g {i — j) ■ It turns out that their proof has a serious gap fixed 
in |KZ] only for the usual Laplacean given by (12.41) . Our presentation, based on an unpublished 
Appendix of [KZj . uses the Laplace method discussed in Subsection 12 . II and is written for admissible 
coupling matrices. Consequently, it works for the hierarchical Laplacean matrix coupling (13. 2p as 
well. 

Reflection Positivity is the missing ingredient. For a basic exposition of the abstract theory see 
[FILSj . If P denotes a plane perpendicular to coordinate axes which divides A into two halves A±: 

A = A+ U A_ and A+ n A_ = 

(P cuts bonds perpendicularly and do not intercept sites of A), let 

r : A+ — ^ A_ 

be a map which assigns to each j G A_|_ its reflected image rj G A_, i.e. the site symmetric with 
respect to P. The reflection map r induces a linear morphism vrp : 2t+ — > 2t_, on the abelian 
algebra Sl-t of bounded function on the configuration space fl± = M^* given by 

By a linear morphism we mean vrp {FG) = Tip {F) iip {G) is satisfied for any F,G & 

Definition 4.3 (Reflection Positivity) A state is said to be a reflection positivity functional 

{F7Tp{F)), > (4.4) 

holds for all F G 21+ . 



According to IFILSI . (■) (n) defined by the expected value with respect to the measure uif^^ 
given by (ILip with A = —A®/, —A the usual Laplacean with periodic boundary conditions (12.41) . 
is a reflection positivity functional. Once (14.41) holds, we have the Schwarz inequality 

(FGnpiFG)), < {Fnp{F)yj'{G7fp{G)yj' . (4.5) 
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Note that the normahzation is unimportant and can be dropped in both sides of the inequahty as 
long as n and are kept fixed. Applying fl4.5l) to every plane P which cuts bonds perpendicularly, 
(■) satisfies the chessboard inequality (see e.g. [FILSj ) 

\{F,{h,)\ <\{l\{F.{h,)\ (4.6) 

\ j I y j \ i I U 

where, for each j, Fj [h, {xi}-^^) = F {h,Xj) is a one parameter h family of bounded function in 
M^. Note that (14. 6 p has a homogenization effect. 

Proof of Theorem \4-^ Repeating the steps of (12. 5p and (12. 6p . we have 

= AT"/^ / n ds, sf-' exp I ^ X: Z'r!"' (1' B) (4-7) 

-^■*+.=i I ^ i=i J 

where 

B = SJS ® / (4.8) 

is the coupling matrix A = J ^ I modified by a matrix S = diag (si, . . . , s„) with the n-vector 
s = (si, . . . , Sn) in the diagonal. Applying the chessboard inequality (14. 6 p to this nonhomogeneous 
partition function (see [KZJ), yields 



zirHi,B)<ll{z!^^ {si A)) 



l/n 



and we have an upper bound 

= (^Vn ^ exv {N hi,'^\s)}J (4.9) 
where, similarly to (12. 7p . 

We are now looking for a lower bound of (14. 7p . Using the O (N) symmetry together with Jensen 
inequality, we have 

» n 

jiN) > TTcis, sf-^Zf)(l,^(J-/iJ)5®/) 

Jlc{l-1/N),c]" -^^ 

= {NRlY" (Zf)(l,5(J-/i/)5®/)> 

> [NRlY'^ Zf) (1, {S ( J - S) ® I) (4.10) 

21 



where (•) denotes the average with respect to the product measure ]^c?/iAr(sj) with d^N^s) = 

ds s^~^xic(i-i/N),c]is) / Rn and c = o^/^a/ (3N^/^ . The constant a is here chosen arbitrarily while 
it has to be tuned properly for the moment generating function. The normalization 



N(3-''I^Rn 

'c(l-l/Af) 

and the two first moments of fi^, 



ds s 



N-l 



1- 1 



N' 



a 1 



1 



(s) = / sd/i7v(s) = V/? 



N1+1/2N 1 _ (1 _ i/AT) 



N+l 



and 



N+l l-(l-l/iV)^ 
iVi+Wi_(l_l/iVf+2 



P 



N + 2 l-(l-l/iV)^ 
are the only quantities that contribute to the lower bound. We have 

{x, {S (J - fil) S)0I x) = {x, {{sf A - (s^) /i) x) 

which, in view of (14.11) . implies 

Zir^l, {S (J - S) ® I) = expQ/iniV(.2)^zf)((.)%A) 

= exp|niv(l;.(.^> + /r)((.n 



Similarly to (12.91) . we have 



^(^) ^ 2»W2-i) (r(iV/2))" 

Ar"(^-i)/ Vdet ( J - /iJ) ■ 



Taking the limit lim \n In^^ /nN , in any order someone wishes, of both fl4.9l) and (14.101) yields 

•n,N^oo 

^f^+llnP + r iP) <-\- ^Eln ( J - /i) < + In + /*(.-2) 



where 



r(/3)= lim 



n,_/V— »oo 

and s is the solution to equation (12.171) . 

Now we solve (12.171) with = /? for /x = and replace this function in the previous equation 
to obtain a lower bound 

rm > fiP) 

with / given by spherical free energy (12.191) . 
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The other side of the same equation gives the following upper bound 
/* < -I In (e/3) + sup ^ (-Eln {J - fi) - Pfi) 

Note that this inequality holds for any /i < 0, so it holds for /i that gives the least upper bound. 
The supremum is attained at = that solves (12.171) with = [3 and we have the upper bound 

/*(/?) < /(/?) 

concluding the proof of Kac-Thompson theorem for admissible reflection positive coupling matrices 
A. 

□ 

The above proof is now modified to establish convergence of the moments generating function. 



4.2 Proof of Theorem 11.11 

The ratio of partition functions 

ef)(Az) = (exp(.(l,=r)))^<., = ^"^li^;"'^'^^ (4.11) 

Zn [p, A) 

where 

Zf)(/3,z,l,A) = y"exp|^(a;,A^) + ^(l,a;)| f[dai,''\x,) (4.12) 

with 1 = ^l/^A^ivj (1, . . . , 1), generates the moments of XnN = (^l/VnN^ We combine 

the procedure of Section fI72\ with the proof of Theorem 14. 21 to evaluate this ratio. Combining (I2.25P 
with (US]) yields 



0nN/2 p n ^ _^ 

/ Y\d''x, exp{—{x,{A-^i)x) + z{l,x] 



: exp 



iV"(^-i)/2 ^det (J-kI) 



2k 



< (v^ T- exp {AT )(.)}) (4.13) 



with 



where the free energy 



/f^ iP,z) = ^\nZr,{P,z;l,A) 
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is a smooth function of {f3,z) such that fi^^ (■5^,0) = fn^'^ (s^) is given by (14.21) . Analogously to 
fl4.10p . we have 

= (iVO"/%xp{nivQ.(.^> + /f) (4.14) 
with the constant a in normalization Rj^ chosen so that 



N 

hm \ -tRn 



2n 



and, recalling equation (I2.15p . the difference between the upper (14.131) and the lower bound (14. 14p 
when the chemical potential k is chosen as a function of f3 is o(l). 

One concludes from the last two equations (I4.13P and (I4.14p the following. If s = s{k,) solves 
the equation (I2.27p . then 

(^^n + hn(3 + (/5, ^)^ < -1 (1 + Eln ( J - k)) - + 0(n) < {^s' + \ns + /f )(.-^ .-)) 

holds provided n and sufficiently large and k = k„ (/?, z) solves 

1 ^2 
(3= E—^ + 



-A — K hNk^ 

for some 0{n) constant independent of z. Substituting exp (^Nf^\l3,z)^ in the numerator and 
denominator (with 2; = 0) of (14. lip , gives 

ei-'tft .) ^ exp {-uN^ f ^ - i!) (1 + „ (1 



and this implies 



lne(/?,z)= lim lne(f)(/?, 



z 



for every /? < /3c, uniformly in compacts of z. 



Remark 4.4 



□ 



1. In order to prove that the theorem holds for the hierarchical Laplacean we need only to prove 
that the functional defined by measure U.l\) satisfies reflection positivity. 
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4.3 Reflection Positivity for the Hierarchical Laplacean 



The 0{N) Heisenberg hierarchical measure has been shown to satisfy reflection positivity by Watan- 
abe who considered the model originally proposed by Dyson with L'^ = 2. We extend his proof to 
the general case L > 1, d > 1. 

Let 6 = {01, . . . , 9k), Ok £ {0, 1}, be the binary representation of a point i G {0, 1, . . . , 2^ — l}. 
With respect to a reflection plane Pk at the k-t\i hierarchy the map r given by 

^"^^ii - \ 1 - if l = k 

acts as exchanging each pair of consecutive blocks of size 2^^"^ indexed by r = (6'a:+i, . . . , 6'x)0 
{ii, . . . ,Z2'=-i, ji, • • • , j2fe-i} — ^ {ji, • • • , j2'=-i,'^i, • • • ,'^2fe-i}- 

For the rf-dimensional lattice, a point i is represented hy = {6i, . . . , 9k) where 9k takes values 
in a box Bl = {0, 1, . . . , L — 1}'^ with periodic boundary conditions: 9k = {9k,i, ■ ■ ■ , 9k,d) with 9k,j 
mod L e {0, 1, . . . , L — 1}. As in the previous subsection, for each hierarchy k a reflection plane 
P is chosen perpendicular to coordinate axes cutting bonds, but not sites, of {0, 1, ... ,L — 1}'' 
dividing this box into two disjoint halves Bl = B^ U B^. With respect to a reflection plane P at 
k-th hierarchy the map r assigns to each 6 its "reflected image" r6 with (r^)^ = 9i if I ^ k and 
(rO), G Bt if 9k G 

Let 

A± = {e: 9k G B^ } 

be the partition of Ak = {O, . . . , — l}'^ into two halves according to a plane P at hierarchy k 
and let ^± denote the set of polynomials in xe, 6 G A±. The reflection map r induces a linear 
morphism vrp : — > given by 



6»GA+ 



To prove reflection positivity of the 0{N) Heisenberg hierarchical measure, its enough to show 
that —A given by (13. 2p is a reflection positivity interaction. For this, let y = B^x G M^^-' and 
notice that y-r G according to whether 9k G Bf^ are coordinate of r = (6';+i, . . . , 9k) for I < k. 
li I > k, then y can be decomposed as 



y = y^ + y , y^ e ^± 



with y = TTpy^. We thus have 

-r:ek&B+ -r.ek£B+ 



Yl ly-l^ + ^'p Y 

T-.e^eBl T:ek(^B+ 



yT? = \\y\\l + Trp WyWl 



^Each pair of consecutive blocks would be reflected : {ii, . . . , 121-1, ji, • . . , j2'=-i} — * {j2'=-ii • ■ • i Jij*2'=-ii • ■ • 1*1} 
if {rd)i — 1 — 61 holds iov I < k. We use the exchange operaction for simplicity. 
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by the morphism property, if / < and 

T 

XX X 

= \\y^\\^ + '^p\\y^\\'^ + '^^yt '^pyt 

T 

and these together with fl3.2p . according to [FILSj . imply that —A is a reflection positivity inter- 
action. 
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